Kinetic theory of shot noise in non-degenerate diffusive conductors 
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We investigate current fluctuations in non-degenerate semiconductors, on length scales inter- 
mediate between the elastic and inelastic mean free paths. We present an exact solution of the 
non-linear kinetic equations in the regime of space-charge limited conduction, without resorting to 
the drift approximation of previous work. By including the effects of a finite voltage and carrier 
density in the contact region, a quantitative agreement is obtained with Monte Carlo simulations by 
Gonzalez et al, for a model of an energy-independent elastic scattering rate. The shot-noise power P 
is suppressed below the Poisson value Ppoisson = 2e7 (at mean current I) by the Coulomb repulsion 
of the carriers. The exact suppression factor is close to 1/3 in a three-dimensional system, in agree- 
ment with the simulations and with the drift approximation. Including an energy dependence of 
the scattering rate has a small effect on the suppression factor for the case of short-range scattering 
by uncharged impurities or quasi-elastic scattering by acoustic phonons. Long-range scattering by 
charged impurities remains an open problem. 

PACS numbers: 72.70.+m, 72.20.Ht, 73.50.Fq, 73.50.Td 
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I. INTRODUCTION 

The kinetic theory of non-equilibrium fluctuations in 
an electron gas was pioneered by Kadomtsev in 1957 
101 and fully developed ten years later The the- 

ory has been comprehensively reviewed by Kogan Q . In 
recent years there has been a revival of the interest in 
this field, because of the discovery of fundamental effects 
on the mesoscopic length scale. (See Ref. jsj for a re- 
cent review.) One of these effects is the sub-Poissonian 
shot noise in degenerate electron gases on length scales 
intermediate between the mean free path £ for elastic 
impurity scattering and the inelastic mean free path 
for electron-phonon or electron-electron scattering. The 
universal one-third suppression of the shot-noise power 
predicted theoretically ||^,^ has been observed in exper- 
iments on semiconductor or metal wires of micrometer 
length |-|ll|. 

The electron density in these experiments is sufficiently 
high that the electron gas is degenerate. The reduction 
of the shot-noise power 



P = 2 



dt' SI{t)SI{t + 1') 



(1.1) 



[with 51{t) the fluctuations of the current around the 
mean current /] below the Poisson value Ppoisson — 2e/ is 
then the result of correlations induced by the Pauli exclu- 
sion principle. When the electron density is reduced, the 
Pauli principle becomes ineffective. One enters then the 
regime of a non-degenerate electron gas, studied recently 
in Monte Carlo simulations by Gonzalez et al. (l2j. In 
a model of energy-independent three-dimensional elastic 
impurity scattering these authors found the very same 
ratio P/Ppoisson = 1/3 as in the degenerate case. The 
origin of the suppression is quite different, however, be- 
ing due to correlations induced by long-range Coulomb 



repulsion — rather than by the Pauli principle. The one- 
third suppression of shot-noise in the computer simula- 
tions required a large voltage and short screening length, 
but was found to be otherwise independent of material 
parameters. 

Subsequent analytical work by one of the authors p^ ] 
explained this universality as a feature of the regime of 
space-charge limited conduction. The kinetic equations 
in this regime are highly non-linear, and could only be 
solved in the approximation that the diffusion term is 
neglected compared to the drift term. This is a question- 
able approximation: The ratio of the two terms is l/d, 
with d the dimensionality of the density of states. The 
resrdt of Ref. pl, 



P/Pp 



12 



3(i2 + 22d + 64 



5 (d + 2)(3d + 4)(3d-f 8)' 



(1.2) 



becomes exact in the large-d limit, when P/Ppoisson 
A/bd, but has an error of unknown magnitude for the 
physically relevant value d = 3. 

The main purpose of the present paper is to report the 
exact solution of the kinetic equations in the space-charge 
limited transport regime. We find that inclusion of the 
diffusion term has a pronounced effect on the spatial de- 
pendence of the electric field, although the ultimate effect 
on the noise power turns out to be relatively small: The 
exact suppression factor differs from Eq. ( |l.2| ) by about 
10% for d = 3. We find 



P/Pi 



p olsson 



0.6857 
0.4440 
0.3097 



for 
for 
for 



d= 1, 
d^2, 
d = 3, 



(1.3) 



close to the values reported by Gonzalez et al. (although 
their surmise that P/ Ppoisson is a simple fraction 1 /d for 
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d = 2, 3 is not borne out by this exact calculation). By 
including the effects of a finite temperature and screen- 
ing length, we obtain excellent agreement with the elec- 
tric field profiles in the simulations (which could not be 
achieved in the drift approximation of Ref. jl^]), and de- 
termine the conditions for space-charge limited conduc- 
tion. We also go beyond previous work by calculating 
to what extent the shot-noise suppression factor varies 
with the energy-dependence of the scattering rate. (This 
breakdown of universality was anticipated in Refs. fl^ ] 
andQ.) 

The paper is organized as follows. The kinetic theory 
is introduced in Sec. ^ where we summarize the basic 
equations and emphasize the differences with the degen- 
erate case. In Sec. ID we formulate the problem for the 
regime of space-charge limited conduction. In Sec. ^ 
we solve the kinetic equations for the case of an energy- 
independent scatte ring rate and compare with the Monte 
Carlo simulations [|2|. We study separately the capac- 
itance fluctuations. The effect of deviations from the 
conditions of space-charge limited conduction is also in- 
vestigated. Energy-dependence in the scattering rate is 
considered in Sec. We conclude in Sec. VI with a dis- 
cussion of the experimental observability in connection 
with electron-phonon scattering. 



The stochastic Langevin current 5J vanishes on aver- 
age, i5J = 0, and has correlator [|| 
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(5J(r, p, i)(5J(r', p', t') = 5{r - r')5{t - t')5{e - e') 



6{h ~ n') / dfi" W,in ■ n")(/ + /" - 2//") 



W,{h-h%f + f'-2ff) 



(2.3) 



determined by the mean occupation number /. (We 
abbreviated /' = /(r,p',t) and analogously for /".) 
The density of states in d dimensions is ^{e) = 
mfl{2meY^^~^ , where we set Planck's constant h = 1. 

A non-degenerate electron gas is characterized by / <C 
1 . In contrast to the degenerate case, the Pauli exclusion 
principle is then of no effect. One consequence is that 
we m ay omit the terms quadratic in / in the correlator 
( ^.3[ ) . A second consequence is that deviations from equi- 
librium are no longer restricted to a narrow energy range 
around the Fermi level, but extend over a broad range of 
s. One can not, therefore, eliminate e as an independent 
variable from the outset, as in the degenerate case. 



B. Diffusion approximation 



II. KINETIC THEORY 

A. Boltzmann-Langevin equation 

Our starting point is the same kinetic theoryR 3 used 
to study shot noise in degenerate conductors ^ ^ . 
We summarize the basic equations, emphasizing the 
differences in the non-degenerate case. The density 
/(r, p, t) of carriers at position r and momentum p — mv 
at time t (where m is the effective mass and v the veloc- 
ity) satisfies the Boltzmann-Langevin equation 



d d . d 

- + v.- + eE(r,i)-^ 



/(r,p,i) + (2.1) 



Here E{r,t) is the electric field (we take the charge of 
the carriers positive), S{r,p,t) is the coUision integral, 
and 6J{r,p,t) is a fluctuating source (or "Langevin cur- 
rent" ) . The collision integral describes the average effect 
of elastic impurity scattering, 

5(r,pn,t) 

= J^W,{h-h')[f{r,ph',t)^f{r,pn,t)]. (2.2) 

The integral over the direction n = p/p of the momentum 
extends over the surface of the unit sphere in d dimen- 
sions, with surface area — 271"^^^ /T{^d). The scattering 
rate W^{n ■ n') depends on the kinetic energy e — /2m 
and on the scattering angle n • n'. The effective mass m 
is assumed to be energy independent. 



We assume that the elastic mean free path is short 
compared to the dimensions of the conductor, so that we 
can make the diffusion approximation. This consists in 
keeping only the two leading terms. 



/(r, nV2me, t) — T{r, e,t) +n- f (r, s, t), 



(2.4) 



of a multipole expansion in the momentum direction n. 
We substitute Eq. (2.4) into the Boltzmann-Langevin 
equation ( [2.l| ) and integrate over fi to obtain the con- 
tinuity equation, 

^p(r, i) + ^ • JCr, £, t) + eE(r, t) ■ ^j(r, e, t) = 0, 

(2.5) 

(2.6) 
(2.7) 



for the energy-resolved charge and current densities 
p{r,e,t) = ei^{e)T{r,e,t), 

j(r,e,t) = -evi^{e){{r,e,t), 
d 



with V — ^J2elm. In the zero- freq uency limit we may 
omit the time derivative in Eq. (2.5). 

Multiplication by n followed by integration gives a sec- 
ond relation between p and j [p2| , 

d 

j(r,e,<) = -D{£)—p{v,£,t) 

- (7(£)E(r, t)^Tir, e, t) + 53{v, e, t), (2.8) 
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a combination of Pick's law and Ohm's law with a 
fluctuating current source. The conductivity a{e) = 
e^v{e)D{e) is the product of the density of states and the 
diffusion constant £>(e) = v^r/d — (2e/md)T(e). The 
scattering rate is given by 



1 



dfi' 



We{-n.-n'){l-n-h'). 



(2.9) 



T(e) J n 
The energy-resolved Langevin current 

JJ(r,e,t) = eT(e)wi/(e) J ^n(5J(r, nV2me, (2.10) 

is correlated as 



SJi{r,e,t)SJ^{r',e',t') 

= 2CT(e)J^(r, £, t)Sim6ir - r')5{t - t')S{e - e'), (2.11) 

where we have omitted terms quadratic in J^. 

These kinetic equations should be solved together with 
the Poisson equation 



d 



(2.12) 



with p{t, t) = J de p{r, e, t) the integrated charge den- 
sity, n the dielectric constant, and peq the mean charge 
density in equilibrium. The Langevin current 5 J induces 
fluctuations in p and hence in E. The need to take the 
fluctuations in the electric field into account selfconsis- 
tently is a severe complication of the problem. 



C. Slab geometry 

We consider the slab geometry of Fig. ||, consisting 
of a semiconductor aligned along the x-axis with uni- 
form cross-sectional area A. A non-fluctuating poten- 
tial difference V is maintained between the metal con- 
tacts at a; = and x = L, with the current source at 
X = 0. The contacts are in equilibrium at tempera- 
ture T. It is convenient to integrate over energy and 
the coordinates r_L perpendicular to the x-axis. We de- 
fine the linear charge density p{x, t) = J dr± p(r, t) and 
the currents I{t) = J dr± J dejx{r,s,t) and SJ{x,t) = 
/ drj^ J de SJx{j^,s,t). The current / is x-independent 
in the ze ro-fr equency limit because of the continuity 
equation (pTsj). We also define the electric field profile 
E{x,t) = J dr^ Ej;{r,t). The vector r_i_ of trans- 
verse coordinates has d — 1 dimensions. The physically 
relevant case is d = 3, but in computer simulations one 
can consider other values of d. For example, in Ref. p^ ] 
the case d = 2 was also studied, corresponding to a hy- 
pothetical "Platland" To compare with the simula- 
tions, we will also consider arbitrary d. 





1 1 ^ 

L ^ 

FIG. 1. Semiconducting slab (grey) between two metal 
contacts (black) at a; = and x = L. The {d— l)-dimensional 
cross-sectional area is A. The current flows from left to right 
in response to a voltage V applied between the contacts. 



For any d the fluctuating Ohm-Fick law (2 
one-dimensional form 



takes the 



Ht) ^--^J J deDie)p{r,e,t) 

+ E{x,t) J dr_L J de J^(r,£,i)^cr(£) 



+ dJ{x,t), 



(2.13) 



where we used that the averages of and E depend on 
X only and neglected te rms quadratic in the fluctuations. 
The Poisson equation ( 2.12 ) becomes 



kA—E{x, t) = p{x, t) - Apcc^, 



(2.14) 



and the correlator (2.11) becomes 



5J{x,t)5J{x',t') 

= 2A5{t-t')5[x-x') I dea[e)T{x,e). 



(2.15) 



Our problem is to compute from Eqs. (2.13) - (2.15) the 
shot-noise power (pTl]). 



D. Energy-independent scattering time 



The Ohm-Fick law (2.13) simplifies in the model of an 
energy-independent scattering time T(e) = r. Then the 
derivative of the conductivity da/de = epv{e) is propor- 
tional to the density of states, and contains the energy- 
independent mobility p — er/m. Eq. (2.13) becomes 



/(O = / "^""^ / de^(e)p(r,£,i) 

+ pp{x,t)E{x,t) +5J{x,t). (2.16) 

The drift term has now the same form ppE as for inelas- 
tic scattering . This simple form does not hold for the 
more general case of energy- dependent elastic scattering. 
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III. SPACE-CHARGE LIMITED CONDUCTION 

For a large voltage drop V between the two metal con- 
tacts and a high carrier density pc in the contacts, the 
charge injected into the semiconductor is much higher 
than the equilibrium charge peq, which can then be ne- 
glected. For sufficiently high V and pc the system enters 
the regime of space-charge limited conduction psf , de- 
fined by the boundary condition 



E{x,t)^0 at 



0. 



(3.1) 



Eq. states that the space charge Q = /g p{x)dx 

in the semiconductor is precisely balanced by the surface 
charge at the current drain. The accuracy of this bound- 
ary condition at finite V and pc is examined in Sec. IVE. 
At the drain we have the absorbing boundary condition 



p{x, t) — at X = L. 



(3.2) 



With this boundary condition we again neglect peq. 

To determine the electric field inside the semiconduc- 
tor we proceed as follows. The potential gain —e(f)(x,t) 
(with E = —d(j)/dx) dominates over the initial ther- 
mal excitation energy of order kT (with Boltzmann's 
constant k) almost throughout the whole semiconduc- 
tor; only close to the current source (in a thin bound- 
ary layer) this is not the case. We can therefore ap- 
proximate the kinetic energy e ~ —e4> and introduce 
this into D{e) and da/de. We assume a power-law 
energy dependence of the scattering time r = roe". 
Then D{e) = (2To/md)e"+i « -{2no/d){-e)°'(j>"+^ 
and dcr/de = (2a + d){To/md)e'^e°'i^{e) w -{2a + 
d){p,o/d){—e)°'^^(l)°'i'{e), wher e we have defined p,o — 
ero/m. Substituting into Eq. ( |2.13|) and using the Pois- 
son equation — kAB^ (j) / dx^ = p we find the third-order, 
non-linear, inhomogeneous differential equation 



(2a- 



dx \dx) dx V dx'^ ) 



2d 



{-eYp.oKA 



[I{t) - 5J{x,t)] 



(3.3) 



for the potential profile 4)(x,t). 

Since the potential difference V between source and 
drain does not fluctuate, we have the two boundary con- 
ditions 



(j){x, t) = at a; = 0, 
<j){x,t) = —V at X — L. 



(3.4) 
(3.5) 



ditions. 



Eqs. (3.1) and ( p.2[ ) imply two additional boundary con- 
{x, t) — at X = 0, 



d_ 

dx 

-^^^(bix.t) — Q at X = L. 
ox^ 



(3.6) 
(3.7) 



We will now solve this boundary value problem for 
(f> = (f) + 5(j), first for the mean and then for the fluctu- 
ations, in both cases neglecting terms quadratic in Scj). 
The case a = of an energy-independent scattering time 
is considered first, in Sec. IV. The more complicated case 
of non-zero a is treated in Sec. 0. 



IV. ENERGY-INDEPENDENT SCATTERING 
TIME 



A. Average profiles 



For a = the averaged equation (O) can be inte- 
grated once to obtain the second-order differential equa- 
tion 



da; 



4-d20 
d dx"^ 



21 



(4.1) 



for the mean potential 4'{x). In this case of an energy- 
independent scattering time r(e) = r we may id entif y po 
with the mobility p = er/m introduced in Sec. II D. No 
integration constant appe ars i n Eq. d4.l| ) , since only then 
the boundary conditions (3.4) and ( |3.6[ ) at x = can be 
fulfilled simultaneously. In Ref. ||l^ the second term on 
the left-hand-side of Eq. (4.1) (the diffusion term) was 
neglected relative to the first term (the drift term) . This 
approximation is rigorously justified only in the formal 
limit d — > 00. It has the drawback of reducing the order of 
the equation by one, so that no longer all boundary con- 
ditions can be fulfilled. Although the soluti on i n Ref. jl^] 
violates the absorbing boundary condition (|3.7D, the final 
result for the shot-noise power turns out to be close to 
the exact result obtained here. 

Before solving this non-linear differential equation ex- 
actly, we discuss two scaling properties that help us along 
the way. Note first that the current / can be scaled away 
by the substitution 



\PkA J 



1/2 



Secondly, each solution x(a;) of 
dx J d da;2 



(4.2) 



(4.3) 



[the rescaled Eq. (4.1)] generates a one-parameter family 
of solutions A"^/2x(x/A). Thus, if we find a solution which 
fulfills the three boundary conditions x(0) — 0, x'(0) — 0, 
x"(l) — (primes denoting differentiation with respect 
to x), then the potential 



/2IL^\ 
\ pnA ) 



1/2 



x{xlL) 



(4.4) 
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solve s Eg . (4.1) with boundary conditions (|3.4l), ( 3.6) 
and (3^). The remaining boundary condition (|3.5| ) de- 
termines the current- voltage characteristic 



^kA ( V 



Vx(l) 



(4.5) 



The quadratic dependence of J on is the Mott-Gurney 
law of space-charge limited conduction . 

We now construct a solution xi^)- One obvious solu- 



tion is Xoix) = aox 



3/2, with 



aa 



1 - 



'Sd 



-1/2 



(4.6) 



This solution satisfies the boundary conditions at x = 0, 
but x'oi^) 7^ for any finite x. Close to the singular 
point X = any solution will approach Xo{x) provided 
that d > 4/3. Let us discuss first this range of d, con- 
taining the physically relevant dimension d — 3. 
We substitute into Eq. (4.3) the ansatz 



1=0 



(4.7) 



consisting of xo{x) times a power series in x^ , with /? 
a positive power to be determined. This ansatz proves 
fruitful since both terms on the left-hand-side of Eq. ( |4.3| ) 
give the same powers of x, starting with order x^ in coin- 
cidence with the right-hand-side. Power matching gives 
Eq. (4.6) for oq and for Z > 1 the conditions 



m=0 



(4- 




0.2 0.4 0.6 0.8 



(4.9) 

The relation with / = 1 is special: It determines the 
power 13, 



(3^+P[2~-d 



'--'-d^O, 
2 8 



(4.10) 



but leaves the coefficient ai as a free parameter [to be 
determined by dema nding that x"(l) = 0]. The positive 
solution of Eq. (4.10) is 



(3^-d-l + -^/9(PT2Id^~32. (4.11) 
8 8 



We find P = (Vl3-l)/4 for d = 2 and /3 = 3/2 for d = 3. 
For / > 2 we solve for ai to obtain the recursion relation 



x/L 

FIG. 2. Profile of the mean electrical potential (j> 
[in units of (27L^/^«:A)^/^], the electric field E [in units 
of (2IL/ij,kA)^^^], and the charge den sity p [in units of 
{2Ik/ ^LAy^^], following from Eq. (|4.l|) for different values 
of d. The drift approximation of Ref. [[f3[ corresponds to the 
case d = 00 in this plot. 



Interestingly enough, the power series terminates for 
d — 12/5, and the solution for this dimension is xi^) — 
-j3/2 _ 13.5/2^ Pqj, arbitrary dimension d > 4/3, the co- 
efficients ai fall off with /, the more rapidly so the larger 
d. We find numerically that the solution with x"(l) = 
has ai = 0.3261 for d — 2 and ai — 0.1166 for d = 3. 



For d < 4/3 we substitute into Eq. (4.3) the ansatz 



ai 



{bu + bio)ao 



(4.12) 



Xix) 



1=0 



■7)/2 



(4.13) 



5 



with 7 = (4 — 3d) /{A — d). Now the coefBcient cq is free. 
Power matching gives further 



Cl 



47(7+1)' 
and the recursion relation 



Cl 



7m - 



^^m — 1 ^IfnCmCl—m 

{dii + d/o)co 



3-7 



— ni) 7m 

a 



(4.14) 



(4.15) 



(4.16) 



for coefficients with I > 2. For d = 1 the solution with 
X"(l) = has Co = 1.3628. 

In Fig. 1^ the profiles of the potential oc the electric 
field E en x' 1 ^^^d the charge density p oc x" are plotted 
for d — 1, 2, 3. We also show the result for d = 00, cor- 
responding to the drift approximation of Ref . |Q . The 
coefficient y(1) appearing in the current-voltage charac- 
teristic (^^) can be read off from this plot. We find 
x(1) = 8/9 for d = 1, x(l) = 0.8180 for d = 2, and 
x(l) = 0.7796 for d = 3. The limiting value for d = cxd is 
X(l) = 2/3. 



The remaining constraint 



-—^ip{x, i) = at X = L 



(4.21) 



(the absorbing boundary condition) will be used later to 
relate 51 to 6J. 

We need the Green function G{x,x'), satisfying for 
each x' the equation C[G{x,x')] = S{x — x'). In view 
of Eq. (|4.3|) for the mean potential one recognizes 



= 3x(a;) - 2x -^x{x) 



(4.22) 



as a solution — which already satisfies Eq. ( 4.19 ). 

Using a standard prescription we find from 'tpi{x) a 
second, independent, homogeneous solution 



-02 (x) = Vi(a;) 



(4.23) 



which fulfills Eq. (4.20). The Wronskian is 



4,^{x)^t^2{x) - i^2{x)^i^i{x) = -x^'^x). (4.24) 
ax ax 

The Green function contains also the factor — 4x/d that 
appears in Eq. ( [1.18 ) in front of the second-order deriva- 
tive of ■0- We find 



B. Fluctuations 

For the fluctuations it is again convenient to work with 
the rescaled mean potential (4.2). We rescale the fluctu- 
ations in the same way, 



{x,t) = - 



21 

IikA 



-1/2 



^p{x,t). 



(4.17) 



We linearize Eq. ( |3.3[ ) with a = around the mean values 
and integrate once to obtain the second-order inhomoge- 
neous linear differential equation 



92 1/, 



da; 



dx^ 

, 5I{t)~5J{x',t) 



dx I dx 



4d^X 
d da;' 



(4.18) 



The integration constant vanishes as a consequence of the 
boundary condition 



0(0;, t) — Q at a; = 



(4.19) 



and the requirement that the fluctuating electric field 
dijj/dx stays finite at a; = 0. [The latter condition actu- 
ally implies d'lp/dx = at a; = 0.] We will solve Eq. ( 4.18| ) 
with the additional condition of a non-fluctuating volt- 
age. 



ipix, t) = at X = L. 



(4.20) 



+ Six' - x)Mx)Mx')] , (4.25) 
where Q{x) = 1 for a; > and 0(x) — for x < 0. 



The solution of the in homo ge neous equation ( |4.1^ 
with boundary conditions (4.19), (4.20) is then 



■ip{x,t) 



dx'G{x,x') / da; 
Jo 



„ 5I{t)~5J{x",t) 



From the extra condition ( 4.21 ) we find 

61{t) =C"^ / dx5J{x,t)g{x), 
Jo 



(4.26) 



(4.27) 



with the definitions 



C-i\(L) £3/2\ , d x'^'jL) 



da; 



xipi (x) 



(4.28) 



, /3x(i) ^\ , d x'/'jL) , Mx') 



(4.29) 



Eq. (4.27) is the relation between the fluctuating total 
current 51 and the Langevin current 5 J that we need to 
compute the shot-noise power. 
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C. Shot-noise power 

Th e sh ot-noise p ower is found by substituting 
Eg. ( 4.27) into Eq. (1.1) and invoking the correlator 
( ^.15 ) for the Langevin current. This results in 



3x(a;)x"(a;) 



X{L) 



ML) 
2L 

(4.41) 



P = 2 / dx 
Jo 



C 



nix), 



nix) ^ 2A / deaie)fix,e). 



(4.30) 
(4.31) 



In order to determine the mean occupation number 
^ix, e) out of equihbrium, it is convenient to change 
variables from kinetic energy e to total energy u = 
s + e(j){x,t). In the new variables x and u we find from 



the kinetic equations (2.5) and (2 
dx 

jix,u) 



-Jix,u) = 
1 



(4.32) 



d 



a[u — ecfiix)]— T ix , u). 
e ax 



The derivatives with respect to x are taken at constant 
u. The solution is 



J^ix,u) = ej(u) 



dx' 



5^(X')] ' 



(4.34) 



where we used the absorbing boundary condition (3.7) 

(which implies .F(L, u) = 0). 

As before [in the derivation of Eq. ( ^3.3| ) from Eq. ( ^.13 )] 
we approximate u — e4>ix) ~ — e0(a;) in the argument of 
CT. (This is justified because < u < kT <^ eV .) Then 
Tix, u) fact orize s into a function of x times a function of 
u, and Eq. ( 4.31 ) gives 



nix)=2elx'-'^ix') dx' x'^-'^ix'), 



(4.35) 



Our final expression for the shot-noise power is 



32 1 d^ 
P = Ppoi.so„^^^y^ dxg^ix)x{x)—x{x). (4.42) 



The scaling properties of x imply that this result does 
not depend on the length L. For d = 1, 2, 3 it evaluates 
to 





f 0.6857 


for 


d = 


1 


(4.33) 


P/Ppoisson = < 0.4440 


for 


d = 


2 




1 0.3097 


for 


d = 


3 



(4.43) 



In Fig. ^ we plot Eq. ( 4.42 ) as a function of the dimension 
d and compare it with the approximate formula ( |1.2| ), 
obta ined in Ref. by neg lectin g the diffusion term in 
Eq. ([4.l|). The ex act r esult ( 4.42 ) is smaller than the ap- 
proximate result ( |1.2| ) by about 10%, 15%, and 25% for 
d = 3, 2, and 1, respectively. For d oo, the drift ap- 
proximation that leads to Eq. ( |l.2[ ) becomes strictly jus- 
tified, and P/Ppoisson approaches 4/5d. The data points 
in Fig. H are the result of the numerical simulation [ p2[ . 
The agreement with the theory presented here is quite 
satisfactory, although our findings do not support the 
conclusion of Ref. |1J] that P = ^Ppoisson in three di- 
mensions. 



where we expressed the result in terms of the rescaled 
potential x- In this equation we recognize the Poissonian 
shot-noise power Ppoisson = 2e/. 



Th e integrals in the expressions (4.28), (4.29), and 
( 4.35 ) for C, g, and n can be performed with th e he lp of 
the fact that x solves the differential equation (4.3). In 
view of this equation. 



d dx \ 



xx' _ 2d 
X'^/^+i ^ ^dd^ 
xipi 4 d 

resulting in 

C = hxiL), 

nix) = Pp. 



(x "''X 



-d/2 /2 



{xx''^ 



XX 



XXX"]X-'"^ 



'oisson^x(a;)x"(a;), 



(4.36) 
(4.37) 
(4.38) 

(4.39) 
(4.40) 




FIG. 3. Shot-noise power P for an energy-independent 
scattering rate as a function of d. The ex act r esult (full curve) 
is compared with the approximate result ( |l2| ) (dashed curve). 
Both curves approach 4/5d for d oo. The data points are 
the results of numerical simulations lll2f . 
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D. Capacitance fluctuations 

The fluctuations 61{t) in the current I{t) are due in 
part to fluctuations in the total charge Q{t) — J dx p{x, t) 
in the semiconductor. The contribution from this source 
to the current fluctuations is 5Iq = {6Q/Q)I. Fluctua- 
tions in the carrier velocities account for the remaining 
current fluctuations Sly = SI — SIq. Since the fluctu- 
ations in Q could be measured capacitatively, it is of 
interest to compute their magnitude separately. Because 
we have assumed that there is no charge present in equi- 
librium in the semiconductor, Q{t) — C{t)V is directly 
proportional to the applied voltage V. The proportion- 
ality constant C{t) is the fluctuating capacitance of the 
semiconductor. (The voltage d oes n ot fluctuate.) 

With the Poisson equation ( 2.14 ) and the boundary 
condition (|^) we have 



The correlator of the capacitance fluctuations, 



Pc = 2 



dtSC{0)dC{t), 



is related to the correlator of 5Iq, 



dt6lQ{0)SlQ{t), 



(4.44) 



(4.45) 



(4.46) 



by Pq = (pIV^ /2kAL)Pc. We also define the correla- 
tors 



Pv = 2 



QV 



dt SIv{0)5Iv{t), 



dt6lQ{0)SIv{t), 



(4.47) 
(4.48) 



such that P = Pq + Pv + Pc 



:qv- 



In view of Eqs. ( |3.3| ), ( 4.1q ) and the boundary con- 
ditions (3.5), (3.7), one obtains E{L) and SE{L,t) as a 



0.01 



o 



o 

Oh 

a 




FIG. 4. Contribution Pq from charge fluctuations to the 
shot-noise power P. The correlator Pc of the capacitance fluc- 
tuations is related to Pq by Pc = {4eK.AL/ jj.V^)PQ /PpoisBon- 

The inte grals can be evaluated by using that x(^) solves 
Eq. (p^ , with the result 



(d + 4)(l-5d) 
2L-3/2 (^) _ I 

a + 4 



(4.56) 
(4.57) 



function of 61 and 6 J, and hence 



In Fig. Q the correlator of the capacitance fluctua- 
tions is plotted as a function of d. For d = 3 we 
find Pc — Q.Q28AeKAL / jiV^ . The corresponding con- 
tribution Pq = 0.0071Ppoisson is relatively small, be- 
ing less than 3% of the contribution from the velocity 
fiuctuations Py — 0.3076Ppoisson- (Incidentally, we find 
that charge and velocity fluctuations are anticorrelated, 
Pqv ~ — 0.0049Ppoisson-) Our calculation thus confirms 
the numerical finding of Ref. jl2j, that the charge fluc- 
tuations are strongly suppressed as a result of Coulomb 
repulsion. However, we do not find the exact cancellation 
of Pq and Pqv surmised in that paper. 



5lM{t) = ]- 1 51 it) - I '^5J{x,t)\ , (4.49) 



5Iv{t) = \I^I{t) + j^'^5J{x,t)^ 



With the help of Eq. ( [4.27D we find 

Pc = i(P + P7-2P,j), 
Pv = \{P + Pj + 2Pij), 

Pcv = \{P-Pj), 
16 



Pij = 



dx{L) 



r dx d^ 

Ppoisson "^^(^)^(^)d^2 

/■^ da; d2 



(4.50) 

(4.51) 
(4.52) 
(4.53) 

(4.54) 
(4.55) 



E. Effects of a finite voltage and carrier density 



For comparison with realistic systems and with com- 
puter simulations one has to account for a finite voltage 
V and a finite carrier density pc in the metal contacts. 
The relevant parameters are the ratios Lc/ L and Lg/ L, 
with Lc — [nkT / epc)^/"^ the Debye screening length in 
the contact and Lg = (kT^/pc)^^^ the screening length 
in the semiconductor. The theory of space-charge lim- 
ited conduction applies to the regime L ^ ^ Lc (or 
kT ^ eV and Pc 3> nV/L"^ — the combination kV/L"^ 
characterizing the mean charge density in the semicon- 
ductor). In this section we will show that, within this 
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regime, the effects of a finite voltage and carrier density 
are restricted to a narrow boundary layer near the cur- 
rent source. We will examine the deviations from the 



boundary condition ( 3A ) and compare with the numeri- 
cal simulations p2[ . 

To investigate the accuracy of the boundary condition 
( |3.l|) we start from the more fundamental condition of 
thermal equilibrium, 



^PcJ^(e)exp(-e/fcT) 
de'u{s')exp{-e'/kT) 



at X = 0. (4.58) 



We keep the absorbing boundary condition p{L, e) = 
at the current drain, since thermally excited carriers in- 
jected from the contact at x = L make only a small 
contribution to the current when eV 3> kT. To simplify 
the problem we assume that all carriers at the current 
source have the same kinetic energy ^dkT, in essen ce re - 
placing the Boltzmann factor exp(— e/fcT) in Eq. ( 4.5g ) 
by a delta function at e = ^kT. We restrict ourselves 
to the physically relevant case d = 3 and subs titut e 
e = ^kT — e(f>{x) in the argument of D{e) in Eq. (2.16). 



Repeating the steps that resulted in Eq. ([f.!]), we arrive 
at the differential equation 



> 
O 



> 
O 




dx 



kT\ d 



21 



e J dx^ iikA 



{x-0- (4.59) 



In comparison to Eq. (4.1), an integration constant ^ ap- 
pears now on the right-hand-side. This constant and the 
current / have to be determined from the four bound- 
ary conditions (f>{0) = 0, K(^"(0) = —pc, 4>{L) = —V, 

4>"iL) = 0. 

We have integrated Eq. ( |4.5E ) numerically. In Fig. || we 
show the electric field for d = 3 and parameters as in the 
simulations of Ref. corresponding to L/Lc — 48.9 
and [Ls/LcY = eV/kT ranging between 40 and 300. We 
find excellent agreement, the better so the larger eV/kT, 
without any adjustable parameter. 



The boundary condition ( |3.1[ ) of zero electric field at 
the current source assumes that the surface charge in the 
current drain is fully screened by the space charge in the 
semiconductor. With increasing eV/kT for fixed L/Lc 
one observes in Fig. ^ a transition from overscreening 
(i? = at a point inside the semiconductor) to under- 
screening {E extrapolates to zero at a point inside the 
metal contact). We c an a pproximate E(x) = — 0g(x — ^), 
where (/jq solves Eq. (4.1) with the boundary conditions 
of space-charge limited conduction. This is an excel- 
lent approximation for eV/kT = 200 {^/L = 0.02) and 
eV/kT = 300 {^/L = -0.004), practically indistinguish- 
able from the curves in Fig. a (top panel). 



FIG. 5. Electric field profiles for eV/kT = 40, 60, 80, 
100, 200, 300, at parameter values d = 3, T = 300 K, 
Pc/e = 10^'' m"^, L = 200 nm, k = 11.7ko (with no the 
diele ctric constant of vacuum) . The full curves follow from 
The data points are the result of numerical simu- 
There are no fitting parameters in this compari- 



Eq. (4.59 



lations 12 



To demonstrate analytically that space-charge limited 
conduction is characterized by the conditions L ^ Ls ^ 
Lc, we will now compute the width of the boundary layer 
and show that it becomes ^ L in this regime. We need 
to distinguish between two length scales ^ and to fully 
characterize the boundary layer. The length ^ determines 
the shift in the asymptotic solution 



0asym(a;) = (l>o{x " C) + 3kT/2e, 



(4.60) 



while the size characterizes the range < x < ^' 
where the exact solution (f>{x) deviates substantially from 

0asyin (*^) ■ 

The values of ^ and ^' are found by comparing Eq. 
(4.60) with the Taylor series 



0(x) 



-Eox^^^+^,^ + 0{x'). (4.61) 

K 2 D 



The coefficients in the Taylor series are determined from 
Eq. (Q, 



° e K pnA' 



(4.62) 
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2 „ kT ^ 21 

6 K e ^kA 



(4.63) 



where 2IJjikA w V'^/L^ up to a coefRcient of order unity 

[cf. Eq. (gj)]. 

We match the two functions ( 4.6C ) and ( 4.61 ) at 
X = S,' , demanding that potential and electric field are 
continuous at a; ~ These two conditions deter- 

mine ^ and Within the regime L Lg ^ Lc we 
find two subregimes, depending on the relative magni- 
tude of Lc/L and {Lg/L)'^. Overscreening occurs when 
LJL > {Lg/L)^. Then Eq « -[2kT en)^/^ , (f)j, « 
(2e/9fcr)i/2(pjK)3/2, and C ~ C = 0{Lc). The differ- 
ence - i = 0{Li/L^) < ^. At the matching point, 
(j) = 0{kT/e), E = OiV^K/pcL^), and p = 0{pc). 
Underscreening occurs when Lc/L <C (i^/i)^. Then 
i?o = 0(l/2«/pL3) « (/.3 = OiplLyn^V'), ^ = 

-0{L^JL^), and C = 0{LI/L^). At the matching point, 
= C'(y4«;3/^3^6)^ £; ^ 0{Eq), and p = 0{pc). In be- 
tween these two subregimes, when L^/L^Lc is of order 
unity, ^' va nishes and ^asym(a;) becomes an exact solu- 
tion of Eq. (4.59) which also fulfills all boundary condi- 



tions. In the same range, ^ changes sign from positive to 
negative values. 

We conclude that the width of the boundary layer 
is of order ma,x{Lc, L^/ L^). At the matching point, 
E V/L. The boundary condition (3T), used to calcu- 
late the shot-noise power P, ignores the boundary layer. 
This is justified because P is a bulk property. We es- 
timate the contribution to -P/Ppoisson coming from the 
boundary layer to be of order max {Lc/L, (Lg/L)'^) (pos- 
sibly to some positive power), hence to be <C 1 in the 
regime of space-charge limited conduction. 



2pqkA 



2/3 



X{x/L). 



(5.1) 



In this way we eliminate the dependence on the current / 
and the length of the conductor L. The rescaled potential 
fulfills the differential equation 



1 -1/2^^ _ 1/2^ ^ ^ 

2 da; dx^ dx^ ' 



(5.2) 



with boundary conditions x(0) = 0, x'(0) — 0, x"(l) = 0- 
We substitute 



(5.3) 



into Eq. (|5.2[). Power matching gives in the first order 
go = 2^^> '^. The second order leaves gi as a free coeffi- 
cient, but fixes the power 77 = (vT3 — l)/2. The coeffi- 
cients gi for I > 2 are then determined recursively as a 
function of gi. From the condition x"(l) = we obtain 
gi = —0.1808. The resulting series expansion converges 
rapidly, with the coefficient gi2 already of order 10"^^. 

The averaged potential and its first and second deriva- 
tive are plotted in Fig. ^. The electric field oc x'{x) 
increases now linearly at the current source, hence the 
charge density oc x"(x) remains finite there. The current- 
voltage characteristic is 



2^oK^ / V 



3/2 



(5.4) 



V. ENERGY-DEPENDENT SCATTERING TIME 

We consider now an energy-dependent scattering time. 
We restrict ourselves to c? = 3 and assume a power- 
law dependence t{s) = tq^". The energy-dependence 
of the rate 1 /r is governed by the product of the scatter- 
ing cross-section and the density of states. For short- 
range impurity scattering the cross-section is energy- 
independent, hence a = —1/2. This applies to un- 
charged impurities in semiconductors. For scattering by 
a Coulomb potential the cross-section is oc e"^, hence 
a = 3/2. This applies to scattering by charged impurities 
in semiconductors The case a = considered so far 
lies between these two extremes We have found an 
exact analytical solution for the case of short-range scat- 
tering, to be presented below. The case of long-range im- 
purity scattering remains an open problem, as discussed 
at the end of this section. 

For short-range impurity scattering, the technical steps 
are similar to those of Sec. IV. We first determine the 
mean potential 4>{x). The scaling properties of Eq. ( |3.3| ) 
are exploited by introducing the rescaled potential xix), 
with 



with x(l) = 0.4559. This is a slower increase of / with V 



than the quadratic increase (4.5) in systems with energy- 
independent scattering. 



0.5 







1 1 

"■■•••p 


1 1 





1 ' 1 







0.2 



0.4 0.6 

x/L 



0.8 



FIG. 6. Profile of the mean electrical potential (p [in 
units of L^(3e^/^//2/io«:A)^''^], the electric field E [in units 
of L{3e^^^I /2fioKA)^^^], and the charge density p [in units of 
{3e^^'^ Ik^^'^ /2poA)^^^] for a three-dimensional conductor with 



short-range impurity scattering, computed from Eq. (5.2) 
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The rescaled fluctuations -0(2;, i), introduced by 



S(f){x, t) 



2/3 



0(x/L,i), (5.5) 



2/ioKyl 

fulfill the linear differential equation 

cm = ij^d^ i^dj^ 

^ 9x3 2 ^1/2 5^2 2 x^/^ dx 
1 ^x'x" , o X'" 



+ 2;t72 



4 y t^3/2 T^i/ 
-(5J(x,t) 



(5.6) 



The solution of the inhomogeneous equation is found 
with help of the three independent solutions of the ho- 
mogeneous equation £[?/'] = 0, 



-02(2;) = x{x) - i^-^x{x), 



VsKx) = -01 (a;) / da; 



(5.7) 



(5.8) 



W2(x') 



where we have defined 

Wix) = i)l(x)'4)'2{x) - '4)'i{x)lp2(x). 



(5.9) 



(5.10) 



The special solution which fulfills '(/'(Oi'^) ~ V''(OjO — 
ip{l,t) = is 



4>ix,t) 



W^ix') 



Six - x')Mx)Mx') 



+ e{x' - x)Mx')Mx) 



V'i(i) 

V'2(l) 



^^„Sm-6_Jix",t)Wix") 



The condition ip"{l,t) — relates the fluctuating cur- 
rent SI to the Langevin current 5J. The resulting ex- 
pression is again of the form (4.27), with now 



C = I dxg{x), 



(5.12) 



^(.) = ^^2 + ^rd.'^^^^(^')^^(-')' 



Xix) \ V2(l) 



W^ix') 



(5.13) 



The shot-noise power is given by Eq. (4.30) with 7i(x) as 
defined in Eq. ( 4.3l|) and the mean occup ation number 
still given by Eq. ( p4| ). Instead of Eq. ( |I35| ) we now 

have 



H(a;) = 2e/x(a;) / dx' 



X{x') 

PpoissonX^^^{x)x"{x), 



(5.14) 



where we integrated with help of Eq. (5.2) and used 

x"(i) = o. 

Collecting results we obtain the shot-noise suppression 
factor 



P/Ppoisson = 0.3777, 



(5.15) 



which is about 20% larger than the result obtained in 
Sec. for an energy - inde pendent scattering time in 
three dimensions. Eq. ( 5.15| ) can be compared with the 
a-dependent result in the drift approximation. 



P/Pi 



Poisson 



6{a - l){a + 2)(16a2 + 36a - 157) 



5(2Q!-5)(8a- 17)(13-f 8a) 



(5.16) 



For a = —1/2 the drift approximation gives P = 
0.4071Ppoisson7 about 10% larger than the exact result 



We now turn briefly to the cas e of long-range impurity 
scattering. The kinetic equation (3.3), on which our anal- 
ysis is based, predicts a logarithmically diverging electric 

1/3 

field (X — In ' a; at the current source for a = 1. In 
the range a > 1, which includes the case a — 3/2 of 
scattering by charged impurities, we could not determine 
the low-x be havi or. [A behavior cx Cx^ is ruled out 
because Eq. (3^) cannot be satisfied with a real coeffi- 
cient C] In the drift approximation, the shot-noise power 
( ^.16 ) vanishes as a — > 1. Presumably, a non-zero answer 
for P would follow for a > 1 if the non-zero thermal en- 
ergy and finite charge density at the current source is 
accounted for. This remains an open problem. 



VI. DISCUSSION 

We have computed the shot-noise power in a non- 
degenerate diffusive semiconductor, in the regime of 
space-charge limited conduction, for two types of elas- 
tic impurity scattering. In three-dimensional systems 
the shot-noise suppression factor P/Ppoisson is close to 
1/3 both for the case of an energy- independent scatter- 
ing rate (P/Ppoisson = 0.3097) and for the case of short- 
range scattering by uncharged impurities (P/Ppoisson = 
0.3777). (The latter case also applies to quasi-elastic 
scattering by acoustic phonons, discussed below.) Our 
results are in good agreement with the numerical simu- 
lations for energy-independent scattering by Gonzalez et 
al. |12]. The results in the drift approximation |13] are 
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about 10% larger. We found that capacitance fluctua- 
tions are strongly suppressed by the long-range Coulomb 
interaction. We discussed the effects of a non-zero ther- 
mal excitation energy and a finite carrier density in the 
current source and determined the regime L ^ Ls ^ Lc 
for space-charge limited conduction (L^ and Lc being 
the screening lengths in the semiconductor and current 
source, respectively). Two subregimes of overscreening 
and underscreening were identified, again in quantitative 
agreement with the numerical simulations [|l2| . 

Let us discuss the conditions for experimental observ- 
ability. We have neglected inelastic scattering events. 
These drive the gas of charge carriers towards local ther- 
mal equilibrium and result in a suppression of the shot 
noise down to thermal noise, P = 8kTdI/dV In- 
elastic scattering by optical phonons can be neglected for 
voltages V < kTjj/e, with Td the Debye temperature. 
Scattering by acoustic phonons is quasi-elastic as long 
as the sound velocity Vs is much smaller than the typ- 
ical electron velocity v « (eV^/m)^/^. For large enough 
temperatures T ^ mvvs/k the elastic scattering time 
r oc e~^/^ depends on energy in the same way as for 
short-range impurity scattering |3l) . 

All requirements appear to be realistic for a semicon- 
ducting sample with a sufficiently low carrier density: 
The electron gas is degenerate even at quite low tem- 
peratures (a few Kelvin). Short-range electron-electron 
scattering is rare due to the diluteness of the carriers. 
Scattering by phonons is predominantly elastic. If the 
dopant (charged impurities) is sufficiently dilute, the im- 
purity scattering is predominantly short-ranged. Under 
these conditions we would expect the shot-noise power to 
be about 1/3 of the Poisson value. 



ACKNOWLEDGMENTS 



We are indebted to J. M. J. van Leeuwen for sho win g 
us how to solve the non-linear differential equation (4.1). 
We thank T. Gonzalez for permission to use the data 
shown in Fig. ||. Discussions with O. M. Bulashenko, 
A. V. Khaetskii, and W. van Saarloos are gratefully 
acknowledged. This work was supported by the Euro- 
pean Community (Program for the Training and Mobil- 
ity of Researchers) and by the Dutch Science Foundation 
NWO/FOM. 



[1] B. B. Kadomtsev, Zh. Eksp. Teor. Fiz. 32, 943 (1957) 

[Sov. Phys. JETP 5, 771 (1957)]. 
[2] Sh. M. Kogan and A. Ya. ShuI'man, Zh. Eksp. Teor. Fiz. 

56, 862 (1969) [Sov. Phys. JETP 29, 467 (1969)]. 
[3] S. V. Gantsevich, V. L. Gurevich, and R. Katilius, Zh. 



Eksp. Teor. Fiz. 57, 503 (1969) [Sov. Phys. JETP 30, 
276 (1970)]. 

[4] Sh. Kogan, Electronic Noise and Fluctuations in Solids 
(Cambridge University, Cambridge, 1996). 

[5] M. J. M. de Jong and C. W. J. Beenakker, in: Mesoscopic 
Electron Transport, edited by L. L. Sohn, L. P. Kouwen- 
hoven, and G. Schon, NATO AST Series E345 (Kluwer, 
Dordrecht, 1997). 

[6] C. W. J. Beenakker and M. Biittiker, Phys. Rev. B 46, 
1889 (1992). 

[7] K. E. Nagaev, Phys. Lett. A 169, 103 (1992). 

[8] F. Liefrink, J. I. Dijkhuis, M. J. M. de Jong, L. W. 

Molenkamp, and H. van Houten, Phys. Rev. B 49, 14066 

(1994). 

[9] A. H. Steinbach, J. M. Martinis, and M. H. Devoret, 
Phys. Rev. Lett. 76, 3806 (1996). 
[10] R. J. Schoelkopf, P. J. Burke, A. A. Kozhevnikov, D. 
E. Prober, and M. J. Rooks, Phys. Rev. Lett. 78, 3370 
(1997). 

[11] M. Henny, S. Oberholzer, C. Strunk, and C. 

Schonenberger, Phys. Rev. B 59, 2871 (1999). 
[12] T. Gonzalez, C. Gonzalez, J. Mateos, D. Pardo, L. Reg- 
giani, O. M. Bulashenko, and J. M. Rubf, Phys. Rev. 
Lett. 80, 2901 (1998); T. Gonzalez, J. Mateos, D. Pardo, 
O. M. Bulashenko, and L. Reggiani, preprint (cond 



mat/9811069) 



[13] 
[14] 
[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 
[22] 

[23] 



[24] 

[25] 
[26] 



C. W. J. Beenakker, pre print (|:ond-mat/9810117D 
K. E. Nagaev, preprint ( |:ond-mat/9812357D . 



M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. B 
51, 16867 (1995); Physica A 230, 219 (1996). 
K. E. Nagaev, Phys. Rev. B 52, 4740 (1995); Phys. Rev. 
B 57, 4628 (1998). 

V. I. Kozub and A. M. Rudin, Phys. Rev. B 52, 7853 
(1995). 

V. L. Gurevich and A. M. Rudin, Phys. Rev. B 53, 10078 
(1995); Pisma Zh. Eksp. Fiz. 62, 13 (1995) [JETP Lett. 
62, 12 (1995)]. 

Y. Naveh, D. V. Averin, and K. K. Likharev, Phys. Rev. 
Lett. 79, 3482 (1997); Phys. Rev. B 58, 15371 (1998); 
Phys. Rev. B 59, 2848 (1999); 

Y. Naveh, Phys. R ev. B 58, R13387 (1998); preprint 
( |cond-mat/9806348| ). 

E. V. Sukhoruko v and D. Loss, Phys Rev. Lett. 80, 4959 
(1998); preprint ( |cond-mat/9809239| ). 
The term rdj/dt is omitted from Eq. ( [2.8[ ) because it 
contributes only to the next higher order in the multi- 
pole expansion. 

A physical realization of the case d = 2 would be a layered 
material in which each layer contains a two-dimensional 
electron gas. A single layer would not suffice because the 
Poisson equation wouid then remai n thr ee-dimensional 
and can not be reduced to the form (2.14). 
In the case of inelastic scattering the drift term fipE fol- 
lows, regardless of the energy dependence of the scatter- 
ing rate, from the equilibrium condition = —kTdJ- / de. 
This condition does not hoid when all scattering is elas- 
tic. 

M. A. Lampert and P. Mark, Current Injection in Solids 
(Academic, New York, 1970). 

N. F. Mott and R. W. Gurncy, Electronic Processes m 



12 



Ionic Crystals (Clarendon, Oxford, 1940). The Mott- 
Gurney law I (xV^ (also known as Child's law) assumes 
local equilibrium and hence requires inelastic scattering. 
Here we find the same quadratic I-V characteristic, but 
only if the elastic scattering time is energy independent 
(cf. Sec. 0). 

[27] E. L. Ince, Ordinary Differential Equations (Dover, New 

York, 1956): § 5.22. 
[28] T. Gonzalez, J. Mateos, D. Pardo, O. M. Bulashenko, 



and L. Reggiani, unpublished data. 
[29] D. Chattopadhyay and H. Queisser, Rev. Mod. Phys. 53, 
745 (1981). 

[30] Formally, one can also consider the case a < —1/2. Na- 
gaev [|l4| has shown that full shot noise, P — Ppoisson, 
follows for a = —3/2. 

[31] V. F. Gantmakher and Y. B. Levinson, Carrier Scatter- 
ing in Metals and Semiconductors (North-Holland, Am- 
sterdam, 1987). 



13 



